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This paper summarizes the yields and the emission patterns of K + and of K~ mesons measured 
in inclusive C+C, Ni+Ni and Au+Au collisions at incident energies from 0.6 A GeV to 2.0 A GeV 
using the Kaon Spectrometer KaoS at GSI. For Ni+Ni collisions at 1.5 and at 1.93 A GeV as well 
as for Au+Au at 1.5 A GeV detailed results of the multiplicities, of the inverse slope parameters of 
the energy distributions and of the anisotropics in the angular emission patterns as a function of the 
collision centrality are presented. When comparing transport-model calculations to the measured K + 
production yields an agreement is only obtained for a soft nuclear equation of state (compression 
modulus Kn ~ 200 MeV). The production of K~ mesons at energies around 1 to 2 A GeV is 
dominated by the strangeness-exchange reaction K~N ^ irY (Y = A, E) which leads to a coupling 
between the K~ and the K + yields. However, both particle species show distinct differences in their 
emission patterns suggesting different freeze-out conditions for K + and for K~ mesons. 



PACS numbers: 25.75.Dw 

I. INTRODUCTION 

Relativistic heavy-ion collisions at incident energies 
ranging from 0.6 to 2.0 A GeV provide a unique opportu- 
nity to study the behavior of nuclear matter at high den- 
sities. These studies are important challenges for testing 
the present understanding of nuclear matter. In addition, 
they are of relevance for astrophysics, as the modelling of 
neutron stars or of supernovae depends on the properties 
of nuclear matter under these extreme conditions [l[ . 

In central Au+Au collisions at the incident energies 
under investigation densities of 2 - 3 times normal nu- 
clear matter density can be reached 0, EJ 0] • A sensitive 
probe to test these conditions is the production of strange 
mesons at or below the production thresholds of these 
particles in free NN collisions. The rest mass of charged 
kaons is 0.454 GeV. For the K + production the threshold 
in NN collisions is 1.58 GeV (in the laboratory system) 
as defined by the associate production NN — > K + AN and 
it is 2.5 GeV for the K~ production via pair creation 
NN — > NNK~K+. 

The key mechanism for the K + production in heavy-ion 



reactions at these incident energies is the accumulation 
of the necessary energy by multiple collisions of particles 
inside the reaction zone. Higher densities increase the 
number of these collisions and especially second genera- 
tion collisions like AN with sufficiently high relative mo- 
mentum to create a K + occur most frequently during the 
high-density phase of the reaction. The density reached 
in the reaction zone depends on the stiffness of nuclear 
matter. Because of their specific production mechanism 
and because of their rather long mean free path (« 5 fm 
at normal nuclear density) K + mesons are ideal probes 
to explore the high-density phase of a heavy-ion reaction 
and to study the stiffness of the nuclear equation of state 

(Eos) iiaii. 

In contrast, the behavior of K _ mesons in a dense nu- 
clear medium is expected to be very different from the 
one of the K + mesons due to two distinct properties: 

(i) The interaction with nuclear matter: The K + are 
hardly absorbed in nuclear matter due to strangeness 
conservation. It is very unlikely that a rare K + (con- 
taining an s quark) encounters an equally rare hyperon 
Y (A, S) containing an s quark. The K+ on the con- 
trary, can easily be absorbed on a nucleon converting it 
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into a hyperon and a pion. Consequently, the mean free 
path of the K~ is significantly shorter than the one of 
the K + . The strangeness-exchange reaction K~N # nY 
has a large cross section and is therefore responsible for 
the appearance and disappearance of K - mesons. This 
channel has been suggested to be the dominant produc- 
tion mechanism in nucleus-nucleus collisions || and this 
has been demonstrated in [l(| [ll|, [ll| • 

(ii) The influence of KN potentials: According to vari- 
ous theoretical approaches the KN interaction is governed 
by the superposition of a scalar and a vector potential 
[lj [H, EDI M, El 01 • While the scalar potential acts 
attractively on both kaon species, the vector potential re- 
pels K + and attracts K~ . For K + these two contributions 
almost cancel leading to a small repulsive K + N interac- 
tion. For the K~ the addition of both attractive interac- 
tions results in a strongly attractive potential. Attempts 
to observe these effects in the res p ective p roduction cross 
sections are under discussion [H [H, Hfl E3, SUM HI • 
These potentials are predicted to have a sizeable effect 
on the azimuthal emission patterns of K + and of K~ (el- 
liptic flow) I25ll26l which has been observed at the Kaon 
SpectrometergJ |2j| . 

This paper intends to give a comprehensive overview 
on the cross sections and on the emission patterns of 
the K + and of the K~ production in mass-symmetric 
heavy-ion reactions in the incident energy range from 
0.6 to 2 A GeV. It summarizes new as well as previ- 
ously published 0, E3, EH H3, Hil results on inclusive 
collisions of Au+Au, Ni+Ni and C+C collisions. Fur- 
thermore, new results focussing on the centrality depen- 
dence of Ni+Ni collisions at 1.5 and at 1.93 A GeV and 
of Au+Au collisions at 1.5 A GeV are presented. Re- 
sults on the kaon production in the mass-asymmetric col- 
lision systems C+Au and Au+C have been published re- 
cently 1 30] . The results on the azimuthal distributions of 
kaons [U |28| as well as of pions and of protons [3l| have 
been published and further publications on this topic are 
in preparation as well as a review on pion production. 

This paper is structured in the following way: First, in 
Section |TT1 a description of the experimental setup and of 
the data analysis is given. In Section Till Al we summarize 
the results on cross sections, energy distributions, and 
polar angle distributions for inclusive collisions, i.e. with- 
out any selection in the collision centrality. Section [III Bl 
presents a detailed study of the centrality dependence of 
the K + and of the K~ production. In Section IIV Al the 
measured yields of K + and of K~ mesons are discussed 
showing that their production is correlated. Despite this 
correlation of the production yields the emission patterns 
of K + and of K~ mesons differ significantly. Section llVBI 
discusses these differences with respect to the influences 
of the KN-potentials and with respect to different emis- 
sion times of the two kaon species. Section llV Cl compares 
the production yields of K + mesons in different collision 
systems to recent transport model calculations to extract 
information on the stiffness of the nuclear equation of 
state. 



II. THE EXPERIMENT 
A. The Setup 

The experiments were performed with the Kaon Spec- 
trometer (KaoS) at the heavy-ion accelerator SIS (Schw- 
erionensynchrotron) at the GSI (Gesellschaft fur Schw- 
erionenforschung) in Darmstadt, Germany. A detailed 
description is presented in [33]. Here we just briefly re- 
view the main features. 

The setup of the quadrupole-dipole spectrometer KaoS 
is shown in Fig. [TJ Positively and negatively charged par- 
ticles are measured separately using different magnetic 
field polarities. The magnetic spectrometer has an ac- 
ceptance in solid angle of f2 sa 30 msr and covers a mo- 
mentum bite of Pmax/Pmin ~ 2. The short distance of 5 
- 6.5 m from the target to the focal plane minimizes the 
number of kaon decays in flight. The loss of kaons by de- 
cay and due to the geometrical acceptance is accounted 
for by corrections which are determined by Monte-Carlo 
simulations using the code GEANT 33]. Particle iden- 
tification is based on momentum and on time-of-flight 
(TOF) measurements. The trigger system is as well 
based on the time-of-flight information to separate pions, 
kaons and protons. For the separation of high momentum 
protons from kaons a threshold Cherenkov detector (34[ 
is used in addition. The trigger system suppresses pions 
and protons by factors of 10 2 and of 10 3 , respectively. 

In total there are three time measurements using seg- 
mented plastic scintillator arrays: The TOF Start detec- 
tor between the quadrupole and the dipole (16 modules), 
the TOF Stop detector in the focal plane of the spec- 
trometer (30 modules), and the Large- Angle Hodoscope 
(LAH) around the target point covering polar laboratory 
angles of 12° < Ub < 48° (84 modules). The latter 
allows for a second time-of-flight measurement for back- 
ground rejection and for the determination of the col- 
lision centrality using the number of measured charged 
particles. 

The trajectory reconstruction is based on three large- 
area Multi-Wire Proportional Counters (MWPC 1 - 
3) [351 ] . one between the quadrupole and the dipole and 
two behind the dipole magnet, each of them measuring 
two spatial coordinates. The efficiencies for kaon detec- 
tion are larger than 95% for each of these detectors. 

The beam intensity is monitored using two scintillator 
telescopes positioned at backward angles (0i ab = ±110°), 
measuring the flux of charged particles produced in the 
target which is proportional to the beam intensity. The 
absolute normalization is obtained in separate measure- 
ments at low beam intensities using a plastic scintillation 
detector directly in the beam line. The beam intensities 
have been chosen such, that the efficiency of the data 
acquisition system (DAQ) due to dead time was always 
above 50%. 

The spectrometer is mounted on a platform which can 
be rotated around the target point on an air cushion in 
a polar angel range from 9\ a b = 0° to 130°. The angular 
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FIG. 1: Top view of the Kaon Spectrometer KaoS with its 
various detector components. 
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interact, prob. 
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3.0 


2.7% 


Ni 


0.8 


2.1% 


Au (0.6,1.0,1-135 A GeV) 


1.0 


3.6% 


Au (0.8,1.5 A GeV) 


0.5 


1.8% 



TABLE I: Thicknesses and interaction probabilities for the 
targets used in the various experiments. 



range covered at each position is A6*i a b = ±4°. Through- 
out this paper we will always quote the mean value. The 
momentum coverage is maximized by measuring different 
magnetic field settings (|-Bdi P oic| = 0.6, 0.9, and 1.4 T). 
The resulting coverage for kaons in rapidity normalized to 
the beam rapidity y/ybcam and in transverse momentum 
p t is sketched in Fig. [5] for three different beam energies 
(1.0, 1.5, and 1.93 A GeV ). The shaded areas correspond 
to different angular settings #i a b of the spectrometer in 
the laboratory as denoted in the figure and to various 
magnetic field settings. 

In this paper we report on measurements of the colli- 
sion systems C+C (0.8, 1.0, 1.2, 1.5, 1.8, and 2.0 A GeV 
beam energy), Ni+Ni (1.1, 1.5, and 1.93 A GeV), and 
Au+Au (0.6, 0.8, 1.0, 1.135, and 1.5 A GeV). The tar- 
gets as well as their respective thicknesses and interac- 
tion probabilities are given in Table U Due to the energy 
loss in the Au targets the average effective beam energies 
-^bcam f° r kaon production in these cases are 0.56, 0.78, 
0.96, 1.1, and 1.48 A GeV. For the other target materials 
the energy loss is negligible. Throughout the text we use 
the values of the nominal beam energies -Ebeam- However, 
in all figures displaying data as a function of the beam 
energy the data points are plotted at E^ am . 

In the case of the 1.5 A GeV Au beam an exceptional 
operation of the GSI accelerator facility was required: ac- 
celeration of 197 Au 63+ ions with the synchrotron SIS up 
to an energy of 0.3 A GeV, then extraction and full strip- 
ping, followed by injection into the Experimental Storage 
Ring (ESR) where the beam was cooled (electron cool- 
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FIG. 2: (Color online) Examples for the coverage of the Kaon 
Spectrometer in transverse momentum p t and in normalized 
rapidity y/ybeam f° r several laboratory angles #i a b (as indi- 
cated) and for various magnetic field settings at 1.0 A GeV, 
at 1.5 A GeV, and at 1.93 A GeV incident energy. 



ing) and finally re-injection into the SIS and acceleration 
up to 1.5 A GeV. 



B. Data Analysis 

1. Track Reconstruction 

For the reconstruction of particle tracks in the spec- 
trometer, reconstruction functions correlating the spatial 
coordinates measured by different detectors are used to 
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combine the information of different detectors to track 
candidates. For example, one reconstruction function 
is used to calculate the x-coordinate in the MWPC be- 
tween the quadrupole and the dipole (x\) as a function of 
the ^-coordinates of the two MWPCs behind the dipole 
(x2 , 2:3) . By comparing a measured position in one de- 
tector with the calculated position based on the hits in 
other detectors using the reconstruction functions, track 
candidates are created. This is done for the spatial coor- 
dinates in the MWPCs as well as for the assignment of 
the modules of the TOF detectors. 

The reconstruction functions are determined by 
Monte-Carlo simulations using a complete description of 
the experimental setup within GEANT. Single tracks are 
followed through the spectrometer and the correlation 
functions are determined by fitting polynomial functions 
up to seventh order and with up to three rr-coordinates 
and up to three y-coordinates to all these simulated 
tracks. 

After the reconstruction a resulting track candidate 
consists of x- and y-coordinates in all the three MWPCs, 
the module numbers, the time and the energy-loss infor- 
mation of the TOF Start and of the TOF Stop detector 
as well as the time and the multiplicity information of 
the Large Angle Hodoscope. 

To determine the efficiency of the tracking procedure 
GEANT is used. This time not only single tracks are gen- 
erated but combinations of one or several tracks and ad- 
ditional background hits in the different detectors as ob- 
served in the experiment. The resulting efficiencies vary 
with the laboratory angle, the magnetic field strength, 
and the size of the collision system because of varying 
track and background multiplicities. The resulting effi- 
ciencies of the tracking procedure are always larger than 
90%. 

For each track candidate constructed in the preceding 
steps, the momentum pi a b and the length of the flight 
path Al between the TOF Start and the TOF Stop de- 
tector modules are obtained from a lookup table gener- 
ated with GEANT. From these quantities together with 
the measured time difference At between the two TOF 
detectors, the squared mass over charge ratio (m/Z) 2 is 
calculated using 
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Since the particles under investigation have Z 
simply use 'mass' in the following. 
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2. Background Reduction and Cross Section Calculation 

Figure [3] shows mass distributions for Ni+Ni collisions 
at 0i a b = 40°. The distribution in the upper panel has 
been measured at a beam energy of 1.1 A GeV and at 
particle momenta pi a b = 0.267 — 0.507 GeV/c using the 
magnetic field setting for positively charged particles. In 



FIG. 3: (Color online) Distributions of the squares of the 
masses assigned to the reconstructed tracks in three differ- 
ent cases: Ni+Ni collisions at 1.1 A GeV (upper panel), at 
1.93 A GeV (middle panel), both using the field polarity for 
positively charged particles, and at 1.93 A GeV for negatively 
charged particles (lower panel). Differently shaded areas show 
the impact of the TOF trigger and of the application of the 
selection criteria during the data analysis. Details see text. 
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this case of low beam energy and of low particle momenta 
the region of the kaon mass is dominated by background. 
The distribution labelled "w/o TOF- Trigger" was mea- 
sured using trigger conditions that forced every particle 
passing through the spectrometer to be recorded. The 
distribution labelled "with TOF- Trigger" shows the clear 
reduction of pions and protons being recorded when us- 
ing the time-of-flight trigger, but still no clear K + signal 
can be seen. The situation is different at higher beam en- 
ergies or at higher particle momenta where the kaons are 
clearly visible already in the distributions "with TOF- 
Trigger" as can be seen in the two lower panels of Fig. [3] 
They show mass distributions for Ni+Ni at 1.93 A GeV 
for positively charged particles in the middle panel and 
for negatively charged particles in the lower panel. In 
the latter case, the trigger condition "w/o TOF-Trigger" 
was not measured. 

To reduce the remaining background effectively with- 
out loosing too many kaons two types of selection criteria 
are applied: 

(i) The so-called "geometrical cuts": These selection 
criteria are based on the comparison between measured 
positions in one of the MWPCs and those extrapolated 
from hits in the other two using the reconstruction func- 
tions described in the previous section. These selection 
criteria are adjusted using measurements during the same 
experiment which are nearly free of background (highest 
beam energy and/or largest laboratory angle). This can 
be done since the particle trajectories inside the spec- 
trometer only depend on the magnetic field and the geo- 
metrical setup but not on quantities like #i a b or -Ebeam- 

(ii) The so-called "velocity cut": For each track can- 
didate the particle velocity between the TOF Start de- 
tector and the TOF Stop detector is calculated as well 
as the velocity between the Large- Angle Hodoscope and 
the TOF Stop detector. The comparison of these two ve- 
locities is a powerful tool to suppress background created 
by fake track candidates. The "velocity cut" is adjusted 
using the measurements nearly free of background as in 
the case of the "geometrical cuts" . 

Depending on the initial signal-to-background ratio the 
strength of these cuts is varied between 5cr and 2a. After 
the application of the selection criteria the final signal-to- 
background ratio varies between 0.7 and 120. The impact 
of the cuts on the mass distributions is as well shown in 
Fig. [3] (labelled "with TOF Trigger and Cuts"). In many 
cases the remaining background is rather small as can be 
clearly seen in the two lower mass distributions measured 
at 1.93 A GeV. Even in the case of low beam energy and 
low particle momenta a clear K + signal can be observed 
after applying the selection criteria (upper panel, please 
note the logarithmic scale). 

In order to subtract the remaining background in the 
mass distributions, a combined fit using a Gaussian and 
a polynomial function to the mass distribution is per- 
formed within a window around the kaon mass. The 
polynomial part is used to estimate the background be- 
low the kaon peak and is subtracted. Figure 0] shows 
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FIG. 4: (Color online) Mass distributions for Ni+Ni collisions 
at 1.1 A GeV at #i a b = 40° for different particle momenta as 
indicated. The solid lines depict the results of a combined 
fit to the background and to the peak at the kaon mass, the 
dashed lines show the background part only. 



mass distributions for Ni+Ni collisions at 1.1 A GeV at 
#iab = 40° as an example with a significant background 
contamination. The fits are performed separately for 
each bin in particle momentum (in most cases having a 
width of 0.05 GeV/c), the upper left graph shows the dis- 
tribution integrated over the full momentum range of this 
particular magnetic field setting for illustration purposes 
only. The solid lines depict the result of the combined fit 
to the background and to the peak at the kaon mass, the 
dashed lines show the background part only. 

The cross sections are calculated from the number of 
kaons iV(p lab , fii ab ) as 



dpi a bdJ?l a b 



iV(pi ab , J?lab) 

1 



target 



/^target ' d 

1 



target " 

N A 

1 



proj face (piab ) ^1. 



ab ) 



e(piab) 



(2) 



with Mt arget being the molar mass, Ptarget being the den- 
sity and dtargot being the thickness of the target material. 
Na denotes Avogadro's constant and N plo j the number 
of projectiles impinging on the target. 

The correction for the geometrical acceptance of the 
spectrometer and for the particle decay / acc (piab, ^lab) is 
calculated using a GEANT simulation. The simulated 
data sets are analyzed with the same analysis procedure 
as the experimental data. The correction is deduced from 
the ratio of particles found after the analysis to those 
initially simulated. 
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FIG. 5: (Color online) Trigger efficiencies etriggor(piab) for 
Ni+Ni at 1.1 A GeV for the three magnetic field settings mea- 
sured (upper panel), and efficiencies of the applied selection 
criteria e cut (piab) (lower panel) for the same data set. 



The total efficiency e(pi a b) is calculated by multiplying 
the detector, the DAQ, and the tracking efficiencies de- 
scribed before as well as the trigger efficiency cogger (Piab) 
and the efficiency of the application of the selection cri- 
teria £ C ut(Plab)- 

To determine the efficiency of the trigger system 
(which is based on particle velocities) pions and pro- 
tons measured at magnetic fields scaled by the mass ra- 
tios TOpion/mkaon or TOproton/^kaon are used which then 
have the same velocities as kaons measured at the nom- 
inal fields. These pions and protons are called "TOF- 
simulated kaons" in the following. Using an open trigger 
condition that forces every event for which a particle en- 
ters the spectrometer to be recorded ("w/o TOF") and 
keeping track of the decisions the TOF trigger would have 
taken, the efficiency of the TOF trigger can be deter- 
mined. The resulting trigger efficiencies are larger than 
90% but they depend on the particle momentum. The 
upper panel of Fig. [5] shows e-trigger(piab) for Ni+Ni at 
1.1 A GeV for the three magnetic field settings measured 
(|-Bdipoie| = 0.6, 0.9, and 1.4 T). 

The efficiency of the applied selection criteria as a func- 
tion of pi a b is determined using the background- free mea- 
surements as described above. In some data sets the 
statistics of these settings was not sufficient. In these 
cases the efficiencies of the "geometrical cuts" were de- 
termined using GEANT simulations, the efficiency of the 
"velocity cut" was determined using the "TOF-simulated 
kaons" . This method has been validated by comparing its 
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FIG. 6: (Color online) Differential cross section as a function 
of piab for K + in Ni+Ni collisions at 1.1 A GeV at #i a b = 
40° . The different symbols denote data measured at different 
magnetic field settings. 



results to efficiencies calculated from real data in cases 
where this is possible. The resulting efficiencies range 
from 75% to 100%. For the example Ni+Ni at 1.1 A GeV 
they are depicted in the lower panel of Fig. [5] 

Figure [H] shows the resulting differential cross sections 
for K + mesons in Ni+Ni collisions at 1.1 A GeV at 
#lab = 40° according to Eq. ((2]) as a function of pi a b- 
This momentum distribution consists of data measured 
at three different magnetic field settings (|-B<iipole| = 0.6, 
0.9, and 1.4 T) which are analyzed separately but overlap 
very well. 



III. EXPERIMENTAL RESULTS 

A. Inclusive Reactions 

In this section the inclusive cross sections for the pro- 
duction of K + and of K _ mesons in C+C, in Ni+Ni and 
in Au+Au collisions are presented as a function of their 
laboratory momentum pi a b as well as of their energy B c .m. 
and their emission angle # c . m . in the center-of-momentum 
frame. Inclusive means that no centrality selection has 
been applied to the data, neither in the analysis nor im- 
plicitly by the experimental setup or the trigger. Let us 
recall that the trigger is generated by a time-of-flight sig- 
nal in the spectrometer and that the multiplicity detector 
is not part of the trigger system. The determination of 
the functional dependencies of the production cross sec- 
tions on E'en, and on # c . m . allows for an extrapolation to 
phase space areas not covered by the experiment and for a 
determination of integrated production yields. All results 
of this section are summarized in Table [TTJ Throughout 
the figures of this paper we use a consistent color code for 
an easy distinction of the collision systems while of course 
keeping full descriptions for black-and-white printed ver- 
sions: Results of Au+Au collisions appear in red, Ni+Ni 
in blue, and C+C in green. 
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Some of the data sets presented in this section are pub- 
lished here for the first time, some have already been 
published previously @, OHO, HH ■ As slightly different 
procedures have been used in these publications to ex- 
trapolate and to integrate the data, we have recalculated 
the total production cross sections using one consistent 
procedure for all data sets. Earlier measurements with 
large errors are not taken into account [36|, 137J, l38[ . 



1. Energy Distributions 

Figure [7] shows the inclusive production cross sections 
for K + (upper panels) and for K~ (lower panels) as a 
function of their momentum pi a b in the laboratory sys- 
tem for three different collision systems. The upper part 
of the figure depicts data taken in Au+Au collisions at 
a beam energy of 1.5 A GeV, the middle part those for 
Ni+Ni at 1.93 A GeV and the lower part those for C+C 
at 1.8 A GeV. To obtain a wide coverage of the phase 
space (see Fig. [2]) measurements at several polar angle 
settings of the spectrometer in the laboratory #i a b have 
been performed. The lines in Fig. [7J are the results of 
simultaneous fits to all angular settings measured for a 
given system assuming a Maxwell-Boltzmann shaped de- 
pendence of the invariant production cross section on the 
center-of-momentum energy E c . m , and a quadratic de- 
pendence on the cosine of the polar emission angle in the 
cm. system c .m.- This procedure will be discussed in 
detail later. 

The necessity for using a non-isotropic distribution in 
6> c . m . is depicted in Fig. [8j Here the invariant cross sec- 

,3 

tions Cinv = are shown as a function of the kinetic 

energy in the cm. system E c . m . — moc 2 for Au+Au at 
1.5 A GeV. Full symbols denote K + , open symbols K~. 
The data measured at a small angle in the laboratory 
system (K+: lab = 32°, K": 6» lab = 40°) are depicted by 
circles, those measured at a large angle (K + : 6^ = 72°, 
K~: #iab = 60°) are represented by squares. The lines 
are Maxwell-Boltzmann distributions 
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chV 
dp 3 



E r 



. exp 



E r 



T 



(3) 



fitted to the data with T being the inverse slope param- 
eter. For an isotropic emission in the cm. system all 
spectra of a given particle type are expected to fall on 
top of each other regardless of the laboratory angle at 
which they have been measured. For K + mesons this 
is clearly not the case pointing towards an anisotropic 
emission. 

Since the distributions have been measured at fixed 
values of #i a b in the laboratory, data points at different 
particle energies correspond to different emission angles 
#c.m. in the cm. system. Therefore, this anisotropy might 
affect the determination of the inverse slope parameter 
T of the energy spectra if data measured at a fixed #i a b 
would be transformed into the cm. system and then fit- 
ted. 
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FIG. 7: (Color online) Inclusive production cross sections for 
K + and for K~ as a function of the laboratory momentum pi a b 
for inclusive reactions of Au+Au at 1.5 A GeV (upper part), 
of Ni+Ni at 1.93 A GeV (middle) and of C+C at 1.8 A GeV 
(lower part). The lines represent a simultaneous fit to all 
laboratory angles using the distribution according to Eq. ([5]). 
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FIG. 8: (Color online) Invariant cross sections for K + (full 
symbols) and for K~ (open symbols) in inclusive Au+Au col- 
lisions at 1.5 A GeV, both at different laboratory angles. The 
lines represent fits according to Eq. <(3j . The observation that 
the data measured at different laboratory angles do not coin- 
cide indicates a non-isotropic polar-angle distribution. 



Therefore, we created "midrapidity distributions" by 
selecting data points within c .m. = 90° ± 10° from mea- 
surements at various laboratory angles. The results are 
shown in Fig. [9l in the upper part for K + , in the lower 
part for K~. The figure summarizes the energy distribu- 
tions of the K + and of the K _ production as measured in 
three different collision systems (Au+Au, Ni+Ni, C+C) 
at different incident energies (0.6 A GeV up to 2 A GeV). 
The measurements of K + in Au+Au at 1.135 A GeV and 
of K~ in C+C at 1.5 A GeV are not displayed in Fig. [5] 
since they do not cover midrapidity. The lines are fits 
to the data according to Eq. ©. The resulting inverse 
slope parameters Tmidrap are given in Table [TTJ 



2. Polar Angle Distributions 

To extract the angular emission pattern we assume 
that the dependence of the invariant cross sections on 
the polar angle c . m . and on the energy -E c . m . can be fac- 
torized. The energy dependence is determined at midra- 
pidity by fitting Maxwell-Boltzmann distributions to the 
data as described above and shown in Fig. [9] As already 
mentioned, each of the data points measured at a given 
center-of-momentum energy E CAri . and a laboratory an- 
gle £?iab corresponds to a different emission angle c .m. in 
the center-of-momentum frame. 

To disentangle the dependencies on the energy and on 
the polar emission angle we normalized each measured 
data point eri nv (-Ec.m., c .m.) to the corresponding value 
owC®c.m.i 0cm. = 90°). The latter is determined us- 
ing the fits to "midrapidity distributions" according to 
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FIG. 9: (Color online) Inclusive invariant cross sections at 
midrapidity as a function of the kinetic energy E c . m . — moc 2 
for K + (upper part) and for K~ (lower part) for the various 
collision systems and beam energies measured. The midra- 
pidity condition is a selection of 9 c .m. = 90° ± 10°. 
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FIG. 10: (Color online) Polar angle distributions for inclusive 
Au+Au collisions at 1.5 A GeV (upper part) and Ni+Ni at 
1.93 A GeV (lower part). Full symbols denote measured data, 
open symbols are reflected at 8 c .m. = 90°. The lines represent 
fits according to Eq. Q. 



ditionally performed simultaneous fits to all momentum 
distributions measured at different laboratory angles 9uh 
for a given system using the combined function 



£-4 = C [1 + a\ 
dp 3 1 1 



cos" 



cm.)] E c 



i. exp 



si 



(5) 

with a 2 , T s f and the normalization C being the three 
variable parameters. The results of this procedure are 
denoted by the super(sub)script 'sf. For Au+Au at 
1.5 A GeV, Ni+Ni at 1.93 A GeV, and C+C at 1.8 A GeV 
the results of these fits are shown as solid lines in Fig. [7] 
The parameters obtained for all collision systems at all 
beam energies are given in Table [TTJ They agree very 
well with the values obtained by the two-step procedure 
denoted by T m id ra p and by a^ lv - The combined fits in ad- 
dition provide the correlations between the three param- 
eters and thus the full error matrix which is necessary for 
calculating the errors of the integrated production cross 
sections. In those cases for which only one angle has 
been measured we take an interval for the polar angle 
anisotropy a| f , denoted by square brackets in Table HT1 
with values set according to the trend at neighboring 
beam energies. This variation of a| f yields additional er- 
rors on the inverse slope parameters T s f as well as on the 
integrated cross sections as tagged by the superscript , &2 , 
in TableE] 



Eq. J3]). Assuming that the energy dependence of the 
kaon production is fully described by these midrapidity 
fits the results are the polar-angle emission pattern. They 
are shown in Fig.[TU]as a function of cos(# c . m .), in the up- 
per part for Au+Au at 1.5 A GeV and in the lower part 
for Ni+Ni at 1.93 A GeV. Full symbols are measured data 
points, open symbols have been reflected at 6* c . m . = 90° 
since for mass-symmetric systems the polar angle distri- 
butions have to be symmetric around # c . m . = 90°. Both 
systems show a forward-backward preference in the emis- 
sion pattern which is more pronounced for K + (left hand 
side) than for K~ (right hand side). 

To quantify the anisotropy, the distributions have been 
fitted with a quadratic dependence on cos(0 c . m .) 

, % , ~ 1 +af cos 2 ((? c . m .) (4) 
dcos(6>c.m.) 

as depicted by the lines in the figure. This procedure has 
been applied to most data sets resulting in the values for 
d2 IV as given in Table [ill In the cases for which only 
one laboratory angle has been measured the coverage in 
0cm. is rather small making the determination of 
impossible. 

The super(sub)scripts 'div' for the angular anisotropy 
and 'midrap' for the inverse slope parameters are used for 
the two-step procedure presented above. Since the energy 
dependence of the kaon production is well described by a 
Maxwell-Boltzmann distribution and the polar angle dis- 
tribution by a quadratic dependence on cos(# c . m .) we ad- 



3. Total Production Cross Sections 

The results of the simultaneous fits have been used to 
extrapolate the data to phase-space regions not covered 
by the experiment and to calculate total production cross 
sections by integrating Eq. ([5]) over the full phase space. 
The extrapolation in E c _ m . contributes with about 35% 
to the total production cross sections. The resulting to- 
tal production cross sections for K + and for K~ for all 
collision systems and for all beam energies are summa- 
rized in Table [TTJ The error bars of the data points in 
the figures showing energy spectra and polar angle dis- 
tributions contain the statistical uncertainties as well as 
point to point systematic errors due to the background 
subtraction. The overall systematic error of the abso- 
lute normalization is stated separately in Table [TTJ and is 
quadratically added to the statistical errors in all figures 
comparing cross section or multiplicities from different 
collision systems. 

In several cases only one polar angle has been mea- 
sured. As already described for the inverse slope param- 
eters T in the previous section, for each of those cases an 
interval for guided by the systematics given by neigh- 
boring beam energies has been used to determine an ad- 
ditional error on the integrated cross section a denoted 
by the superscript 'a2' in Table ITTI 

Part of the data has been published earlier with the 
methods to extrapolate the measured data to the full 
phase space being slightly different in the various publi- 
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cations. In this paper, we consistently apply one single 
method to extrapolate and integrate the data. 

For C+C collisions a and T for K + and for K~ were 
published [19] assuming the angular anisotropy for all in- 
cident energies to contribute with 20% to the total cross 
section independent of the incident energy which is equiv- 
alent to a,2 = 0.6. Now we determine a 2 for each mea- 
surement separately. The differences between the values 
for a and T in Table [IT] and those published in [lj| are 
nevertheless smaller then the statistical errors. 

In the case of Ni+Ni, cross sections for the K + and for 
the K~ production at 1.93 A GeV have been published 
in [20I ] with the angular anisotropy being taken from a 
two-step procedure rather than from a simultaneous fit. 
Also in this case the differences between the results in 
Table [IT] and the previously published values are smaller 
than the statistical errors. 

For Au+Au a, T and 02 have been published for 
K+ 0]. The results in Table [n] for 0.6, 0.8, 1.0 and 
1.135 A GeV differ less than the statistical errors from 
the published values. For K + at 1.5 A GeV the results 
published in Q are: a = 267 ± 30 mb, T = 100 ± 5 MeV 
and 02 = 1.06 ± 0.3. They correspond to a low statistics 
measurement using a thicker target reducing the effective 
beam energy for the K + production to 1.46 A GeV. This 
difference in effective energy accounts for a difference of 
about 10% in cross section compared to the new high 
statistics experiment reported upon in this paper using 
a thinner target (with E e ff = 1.48 A GeV). In addition 
in the experiment described in [jj the anisotropy in the 
polar angle emission pattern was underestimated due to 
a reduced coverage in 9\ a h accounting for an additional 
5% difference in the cross section. 

To calculate particle multiplicities M the integrated 
production cross sections a (see Table Hi]) need to be di- 
vided by the total reaction cross section ay which can- 
not be determined easily as the particle-multiplicity dis- 
tribution measured with a multiplicity trigger condition 
has a cutoff at low multiplicities (see also Sect. IIIIBI) . 
Therefore, we determine the reaction cross sections us- 
ing Glauber calculations [3!| resulting in ay (Au + Au) = 
6.8 barn, a r (Ni + Ni) = 3.1 barn, and ay(C + C) = 
0.95 barn. Figure [TT1 summarizes the multiplicities of K + 
and of K~ mesons as a function of the beam energy as 
determined in inclusive reactions (Au+Au, Ni+Ni, C+C) 
and normalized to the mass number A of the respective 
colliding nuclei. Both particle species exhibit strongly 
rising excitation functions as expected due to the prox- 
imity of the thresholds in binary NN collisions (1.58 GeV 
for K+, 2.5 GeV for K"). The solid lines reflect fits ac- 
cording to the formula [40j 



Mk 
A 



exp 



tin- 



(6) 



with Tmax = To • (-Ebcam)' 7 - The variable parameters in 
the fit are C, T and ij. As well for K + as for K~ the 
multiplicities per mass number A increase with system 
size from C+C to Au+Au at the same incident energy. 



10 



10 



10 



K + K 




beam 



[AGeV] 



FIG. 11: (Color online) Multiplicities of K + (full symbols) 
and of K~ mesons (open symbols) per mass number A of the 
respective collision system as a function of the beam energy. 
The lines represent fits to the data according to Eq. ([6]). 



B. Centrality Dependence 



The collision centrality was derived from the multi- 
plicity of charged particles measured in the Large- Angle 
Hodoscope (Mullah)- Figure fT2l shows the respective 
multiplicity distributions for Au+Au at 1.5 A GeV (up- 
per panel) and for Ni+Ni at 1.93 A GeV (lower panel) 
measured with a multiplicity trigger. In order to study 
the centrality dependence of the K + and of the K~ pro- 
duction the data were grouped into five centrality bins 
both for Ni+Ni at 1.5 and at 1.93 A GeV as well as 
for Au+Au at 1.5 A GeV. These bins are also depicted 
in Fig. 1121 The distributions have been normalized to 
the beam intensity, to the target thickness and to the 
efficiency of the DAQ system so that the area between 
the respective bin boundaries represents the correspond- 
ing fraction of the total reaction cross sections o r for 
a given bin. Very peripheral collisions in the first cen- 
trality bin might be missed by the multiplicity trigger, 
however, kaons from such peripheral events are measured 
since they are triggered by the time-of-flight trigger in 
the spectrometer. The fraction of a r for this most pe- 
ripheral bin is determined by taking the inclusive total 
reaction cross section from a Glauber calculation [3^| 
and subtracting the sum of the experimentally measured 
values of the four other centrality bins. For Au+Au, 
the five bins correspond to - 5.4%, 5.4% - 18.1%, 
18.1% - 31.1%, 31.1% - 52.3%, and 52.3% - 100% 
of a r from central to peripheral collisions, for Ni+Ni to 
- 4.4%, 4.4% - 15.0%, 15.0% - 26.5%, 26.5% - 45.9%, 
and 45.9% — 100.0% of a r . The corresponding mean 
numbers of participating nucleons j4 par t have as well been 
calculated using Glauber calculations. 
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reaction 



rani 

[A GeV] 



a 
[mb] 



[MeV] 



sf 

a ■, 



midrap ( 3'2 

[MeV] 



K+ 



c+c 


0.8 


44 


0.016 ± 0.002 stat ± 0.0016 syst ± 0.002 a2 


54 ± 4 ± l a2 


[-0.5, 0.5] 


52 ±4 




1.0 


44,54,70 


0.1±0.01 stat ±0.01 syst 


60 ±3 


0.25 ±0.27 


56 ±4 




1.2 


40 


0.3 ± 0.03 stat ± 0.03 syst ± 0.02 a2 


67 ± 4 ± 2 a2 


[0.0, 1.0] 


67 ±5 




1.5 


32,48 


1.3±0.1 stat ±0.13 syst 


77 ±3 


0.67 ±0.30 


77 ±5 




1.8 


32,40,48,60 


3.15±0.13 stat ±0.32 syst 


81 ±2 


1.21 ± 0.15 


81 ±3 




2.0 


32,40 


5.1 ± 0.3 3tat ± 0.5 syst 


86 ±3 


1.13 ± 0.23 


85 ±5 


Ni+Ni 


1.1 


40 


5.9 ± 0.3 stat ± 0.6 syst ± 0.3 a2 


87 ± 3 ± l a2 


[0.5, 1.0] 


87 ±4 



1.5 40 31.4 ± 1.3 stat ±3.1 syst ± 1.7 a2 101 ± 3 ± 3 a 

1.93 32,40,50,60 95 ± 2 stat ± 9.5 syst 112 ±2 



[0.5,1.0] 97 ±7 
0.90 ± 0.06 108 ± 5 0.85 ± 0.06 



Au±Au 



K" 

C±C 



0.6 
0.8 
1.0 
1.135 
1.5 



50 
44,84 
44,84 

56 

32,40,48,60,72 



0.5 ± 0.1 stat ± 0.05 syst ± 0.03 a 
8.7 ± 1.4 stat ±0.9 syst 
31 4 ± 4 stat ± 3 \ a y at 



65 ± 14 s 



346 ± 9 s 



6.5 syst ± 2 a2 
± 35 syst 



49 ± 6 ± l a2 
67 ±4 
82 ±3 

89 ± 9 ± l a2 
111 ± 2 



[0.75, 1.25] 
1.20 ± 0.42 
1.07 ±0.24 
[0.75, 1.25] 
1.25 ± 0.09 



50 ±7 

66 ±5 1.25 ±0.4 
81 ±6 1.16 ±0.2 

111 ±5 1.28 ±0.09 



1.5 40 0.016 ± 0.006 stat ± 0.0016 syst ± 0.001 a2 

1.8 40,60 0.078 ± 0.01 stat ± 0.008 syst 

2.0 40 0.15 ± 0.03 stat ± 0.015 syst ± 0.004 a2 



50±12±7 a2 [0.0,1.0] 

65 ±8 0.80 ± 0.25 55 ± 8 
57±7±2 a2 [0.5,1.5] 58 ±9 



Ni±Ni 


1.5 


40 


0.58 ± 0.05 stat ± 0.06 syst ± 0.03 a2 


89 ± 6 ± 2 a2 


[0.5, 1.0] 


98 ± 10 






1.93 


32,40,50,60 


2.6±0.1 stat ±0.26 syBt 


89 ±2 


0.70 ± 0.09 


84 ±4 


0.66 ±0.09 


Au±Au 


1.5 


32,40,48,60 


5.0 ± 0.4 stat ± 0.5 syst 


87 ±4 


0.56 ±0.26 


91 ±8 


0.64 ±0.26 



TABLE II: Summary of the integrated production cross sections of K + and of K~ mesons. The inverse slope parameters T B { 
and Tmidrap as well as the angular anisotropics a| and a^ lv have been determined by two different procedures as explained 
in the text. For systems with only one measured angle #i a b the angular anisotropics a| f have not been fitted but set to an 
interval guided by neighbouring beam energies. It is given in square brackets. The additional errors for a and T s j m due to this 
procedure are denoted by the superscript 'a2'. 



1. Multiplicities 

The multiplicity of a particle species for each centrality 
bin is denned as M = <r/(/ • cr r ) with a being the pro- 
duction cross section for the respective particle species 
and (/ • <7 r ) being the fraction of the total reaction cross- 
section for the respective bin. Figure fTSI presents the mul- 
tiplicities per number of participating nucleons M/A paTt 
as a function of A palt at a beam energy of 1.5 A GeV for 
Au±Au and for Ni+Ni. To calculate these multiplicities 
only data measured at 6>i a b = 40° have been used since 
for Ni+Ni this is the only laboratory angle measured at 
that beam energy. The K + are shown in the upper panel, 
the K~ in the middle panel and the pions in the lower 
panel. Since neutral particles can not be detected by the 
spectrometer, the total pion multiplicity has been calcu- 
lated as M(tt) = 3/2M(tt+) + 3/2M(tt~). In all panels 
full symbols denote Au+Au, open symbols Ni+Ni. 

The lines in Fig. [13] are functions M ~ ^4p ar t fitted 
to the data separately for Au+Au (continuous lines) and 
for Ni+Ni (dashed lines). For both systems the pions 
show a linear dependence of M on A part with a close to 1 



(a^Au) = 0.96 + 0.05, a„ (Ni) = 1.0 + 0.05) which means 
that the number of pions produced is proportional to the 
number of nucleons participating in the reaction. For K + 
on the other hand the multiplicities rise stronger then 
linear with A part (a K +(Au) = 1.34 + 0.16, a K +(Ni) = 
1.26 ± 0.06). The same holds for K" (a K - (Au) = 1.22 ± 
0.27, a K -(Ni) = 1.25 + 0.12). 

Both, the K + and the K~ multiplicities rise similarly 
with centrality leading to a nearly constant K~ /K + ratio 
as a function of A palt (see Fig.[l8]in Sect. lIV A| although 
their production thresholds differ significantly. This will 
be discussed in detail in Sect. IIV Al 

2. Energy Distributions 

Although the multiplicities of K + and of K~ mesons 
per j4 par t in Au+Au and Ni+Ni show the same rise with 
^part and even have the same absolute values, significant 
differences between K + and K~ have been found [lol ]. 

Figure [14] shows the invariant cross sections for K + 
and for K - mesons measured close to midrapidity (#i a b = 
40°) as a function of the kinetic energy in the center-of- 
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FIG. 12: The reaction cross section as a function of the 
charged particle multiplicity in the Large-Angle Hodoscope 
for Au+Au collisions at 1.5 A GeV (upper panel) and for 
Ni+Ni at 1.93 A GeV (lower panel). The shaded areas denote 
the five centrality bins. 



momentum system, in the upper part for Au+Au col- 
lisions at 1.5 A GeV, in the lower part for Ni+Ni at 
1.93 A GeV. The left panel of each graph depicts K + , 
the right panel K~. The uppermost distributions corre- 
spond to the most central reactions, the subsequent bins 
are shown from the top to the bottom of each graph with 
decreasing centrality. The solid lines represent Maxwell- 
Boltzmann distributions according to Eq. fitted to 
the data. 

The resulting inverse slope parameters T for K + and 
K~ mesons are shown in Fig. [15] as a function of 
^4part- The upper part displays Au+Au and Ni+Ni col- 
lisions at 1.5 A GeV, the lower part Ni+Ni collisions at 
1.93 A GeV. In Fig. [T5l two distinct features can be ob- 
served: 

(i) The inverse slope parameters increase with central- 
ity for both particle species for both collision systems and 
both beam energies. 

(ii) The inverse slope parameters of the K + spectra are 
larger than those of the K~ spectra for both systems, at 
both energies and for all centralities. This is discussed in 
detail in Sect. ITVBl 



3. Polar Angle Distributions 

Another observable showing a distinct difference be- 
tween K + and K _ is their polar angle emission pattern. 
Deviations of the angular distributions from isotropy 
can be studied by the ratio <7i n v(6> c . m .)/<7i n v(90 o ) as a 




100 200 300 
A 



part 



FIG. 13: (Color online) Dependence of the multiplicities of 
K + (upper panel) and of K~ mesons (middle panel) as well as 
of pions (lower panel) on A part . Full symbols denote Au+Au, 
open symbols Ni+Ni, both at 1.5 A GeV. The lines are func- 
tions M ~ j4p ar t fitted to the data separately for Au+Au 
(solid lines) and Ni+Ni (dashed lines). The data have been 
measured at #i a b = 40°. 



function of cos(# c . m .) as demonstrated in Sect. Ill A. 
Here, crinv(0c.m.) are the invariant particle production 
cross sections measured at polar angles 9 c .m. in the 
center-of-momcntum frame and <7i nv (90°) is deduced 
from Maxwell-Boltzmann fits to the midrapidity spec- 
tra. These spectra have been obtained for the centrality- 
dependent data in the same way as described for the in- 
clusive data in Sect. IIII Al For an isotropic distribution 
this ratio would be constant and identical to 1. 

Due to limited statistics we divided the data sets into 
two centrality bins only: non-central collisions (15 — 
100 % of cr r for Ni+Ni, 18.1 - 100 % of cr r for Au+Au) 
and near-central collisions (0 — 15 % of <7 r for Ni+Ni, 
- 18.1 % of cr r for Au+Au). 

Figure [16] displays the polar angle distributions for 
Au+Au at 1.5 A GeV (upper part) and for Ni+Ni at 
1.93 A GeV (lower part). The upper panels of both parts 
of the figure show K + , the lower panels K~, both for non- 
central (left) and near-central collisions (right). 

As for the inclusive data the distributions have been 
fitted with a quadratic dependence on cos(# c . m .) accord- 
ing to Eq. ^ to quantify the anisotropy. The solid lines 
and the values for the parameter a 2 in Fig. \W\ represent 
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FIG. 14: (Color online) Energy distributions of the invariant 
cross sections for Au+Au at 1.5 A GeV (upper part) and for 
Ni+Ni at 1.93 A GeV (lower part) close to midrapidity for dif- 
ferent centralities. The uppermost spectra correspond to the 
most central collisions. The subsequent bins are shown from 
top to bottom with decreasing centrality. The lines represent 
fits to the data according to Eq. ([3]). 
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FIG. 15: (Color online) Inverse slope parameters of the energy 
distributions of K + and of K~ at 1.5 A GeV (upper part) and 
at 1.93 A GeV (lower part) as a function of A part . 



IV. DISCUSSION 

In this section we discuss the observed centrality and 
system-size dependences of the K + and of the K~ multi- 
plicities in the context of the production mechanisms of 
these particles, we address the dynamics of the emission 
of the two particle species in terms of energy spectra and 
polar angle distributions, and we extract information on 
the stiffness of the nuclear equation of state by compar- 
ing the K + multiplicities from different collision systems 
to recent transport model calculations. 



the results of these fits. For both systems the K~ mesons 
exhibit a nearly isotropic emission pattern in near cen- 
tral collisions whereas the emission of the K + mesons is 
forward-backward peaked. 



A. The Connection between the K and the K + 
Production 

As presented in Sect. IIIIBl the production yields of K + 
and of K~ mesons exhibit a very similar dependence on 
the collision centrality. Figure [13] shows that the multi- 
plicities of both kaon species exhibit the same rise with 
the number of participating nucleons ^4 par t despite the 
fact that the thresholds for the production of the two 
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FIG. 16: (Color online) Polar angle distributions as mea- 
sured in Au+Au at 1.5 A GeV (upper part) and in Ni+Ni 
at 1.93 A GeV (lower part). Full symbols denote measured 
data, open symbols are reflected at # c .m. = 90°. The lines 
represent fits according to Eq. ((4}. 



particles species in binary NN-collisions are very differ- 
ent. This is observed in Au+Au as well as in Ni+Ni 
collisions. 

Figure [T7] shows the multiplicities of K + mesons from 
inclusive reactions as a function of the system size A at 
several incident energies as well as those of K~ mesons at 
1.5 A GeV. To interpolate between measured data points 
in case of slight differences in the effective beam energies 
due to different energy losses in the respective targets 
the fits to the excitation functions according to Eq. (|6|) 
as shown in Fig. [TT] have been used. The lines in Fig. [T7] 
are functions M ~ A 1 fitted to the data with the re- 
sulting values for 7 given in the figure. Please note: To 
distinguish the two approaches, we use the exponent 7 to 
quantify the rise of the multiplicities from inclusive reac- 
tions as a function of the system size A and the exponent 
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FIG. 17: (Color online) Multiplicities per mass number of the 
collision system M/A as a function of A for C+C, Ni+Ni, and 
Au+Au. The lines represent the function M ~ A 1 fitted to 
the data. 



a for the rise with the number of participating nucleons 
Apart as determined from the analysis of the centrality 
dependence in Sect. IIII Bl 

For K+ at 1.5 A GeV 7 = 1.32 ± 0.06 is extracted 
which is almost identical to the value of a = 1.34 ± 0.16 
as determined from the dependence of the multiplicity M 
on Ap ar t for the K + production in Au+Au at the same 
energy as shown in Fig. [13] As can be seen in Fig. [17] 7 
increases towards lower incident energies. This is in good 
agreement with the assumption of multiple collisions be- 
ing needed to accumulate the necessary energy for the 
K + production at beam energies below the threshold in 
binary nucleon-nucleon collisions. The lower the incident 
energy the more collisions seem to be needed. Since the 
densities reached in heavy reaction systems are signifi- 
cantly larger than in light systems the difference between 
the K + yield in C+C and Au+Au increases with decreas- 
ing beam energy. The same holds for the differences be- 
tween peripheral and central collisions in which different 
levels of baryon densities are created. The similarity be- 
tween the dependencies on A part and on A is therefore not 
astonishing. These observations will be used in Sect. II V Cl 
to extract the stiffness of the nuclear equation of state. 

The yield of the K~ mesons as a function of the system 
size A increases with 7 = 1.31 + 0.11 at 1.5 A GeV . This 
is roughly the same rise as for the dependence on A part 
in Au+Au which was determined to be a — 1.22 ± 0.27 
(see Sect. IIII Bl) . As in the case of the A part dependence 
the rises of the K~ and of the K + multiplicities with 
A are rather similar although the K production hap- 
pens much farther below its respective nucleon-nucleon 
threshold (2.5 GeV) than the K+ production (1.58 GeV). 
In the case of the analysis of the centrality dependence 
at 1.5 A GeV not only the same rise of the K + and of the 
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FIG. 18: (Color online) K /K + ratio as a function of Apart 
for Au+Au at 1.5 A GeV (full circles) and for Ni+Ni at 
1.5 A GeV (open squares) as well as at 1.93 A GeV (full 
squares). The dashed lines denote the K~/K + ratios as cal- 
culated within the statistical model [4lj . 
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collisions as a function of the beam energy. 



K multiplicities with centrality was observed but also a 
rough agreement between the data measured in Au+Au 
and Ni+Ni (see Fig. [13]). 

The similar rise with A part yields a rather constant 
K~/K + ratio as a function of A palt as can be seen in 
Fig. [TH] for three cases. At 1.5 A GeV the ratios for 
Au+Au and for Ni+Ni are the same. The dashed lines 
denote the K~/K + ratios as calculated within the sta- 
tistical model [41] . 

The similar rise of both, K + and K~, as a function of 
the collision centrality as shown in Fig. [13] and as a func- 
tion of the system size (see Fig. [T7]) suggests that the 
production mechanisms of the two kaon species might 
be linked. As already suggested i n H and supported by 
transport model calculations [H, the K in heavy- 
ion collisions at SIS energies are mainly produced via 
the strangeness-exchange reaction 7rY # K~N, with Y 
denoting the hyperons A and E. On the other hand, 
strangeness has to be conserved when producing these 
hyperons and the energetically most favorable way is to 
produce them together with K+ (and K°) mesons (asso- 
ciate production). Thus the production of K + and of K~ 
mesons is coupled via the strangeness-exchange reaction 
and the K - inherit the same dependence on the system 
size and on the collision centrality. In [12] it was argued 
that the strangeness-exchange channel reaches chemical 
equilibrium resulting in the K~/K + ratio to be propor- 
tional to the pion density and that such a proportionality 
was observed for beam energies lower than approximately 
10 A GeV as reached at the SIS and at the AGS (Al- 
ternating Gradient Synchrotron) accelerators at the GSI 
and at the Brookhaven National Laboratory (BNL). 



B. The Dynamics of the K + and of the K 
Emission 

The strangeness exchange reaction couples the yields of 
the K~ and of the K + mesons as discussed in Sect. II V A"l 
On the other hand, K + and K~ show rather distinct dif- 
ferences in observables like energy spectra or polar angle 
distributions which are sensitive to the dynamics of the 
particle emission. 

FiguTeHSlpresents the inverse slope parameters T s f (see 
Table [IT]) as determined by the simultaneous fits to the 
momentum distributions of the production cross sections. 
They are shown as a function of the incident energy for 
inclusive C+C, Ni+Ni and Au+Au reactions, on the left 
hand side for K + , on the right hand side for K~. The 
inverse slope parameters are higher for heavier systems 
and in case of the K+ they rise with increasing incident 
energies. They are always higher for K + than for K~ 
mesons at the same beam energy. 

The same trend is as well observed as a function of 
the collision centrality as shown in Fig. [15] in Sect. IIII Bl 
Figure [2UJ shows the correlation between the inverse slope 
parameters T(K~) of the K~ mesons and those of the 
K + mesons T(K + ) measured in the same collision system 
and at the same incident energy. For C+C the results 
from inclusive collisions at 1.5, 1.8, and 2.0 A GeV are 
shown. For Ni+Ni at 1.5 and at 1.93 A GeV as well 
as for Au+Au at 1.5 A GeV the results of the centrality- 
dependcnt analyses are shown. The dashed line indicates 
equal inverse slope parameters for K - and for K + . The 
data clearly deviate from this line. The inverse slope 
parameters of the K + are about 15 - 25 MeV higher than 
those of the K~, independent of the collision centrality, 
of the collision system and of the beam energy (within 
the measured energy range from 1.5 to 2.0 A GeV). 
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FIG. 20: (Color online) Relation between the inverse slope 
parameters of the energy distributions of K~ and of K + in 
various collision systems and at different incident energies (see 
text for details). 



While the measured K~/K + ratio agrees well with 
statistical- model calculations [U] , the different values for 
the inverse slope parameters for K + and for K~ clearly 
contradict the assumption of a simultaneous freeze out 
of both kaon species. For K~ mesons the chemical and 
the kinetic freeze out coincide as nearly no elastic scat- 
tering occurs due to the strong absorption. If the kinetic 
decoupling of the K + is at a higher "temperature" as the 
chemical freeze out of the K - , they cannot have a unique 
chemical decoupling. 

Different inverse slope parameters for K + and for K - 
mesons have been proposed to result from the influ- 
ence of the repulsive and of the attractive KN poten- 
tials in early Boltzmann-Uehling-Uhlenbeck (BUU) and 
in Relativistic-Quantum-Molecular-Dynamics (RQMD) 
transport-model calculations [U, [H, |44 |. The experi- 
mentally observed difference between the inverse slope 
parameters of K + and of K~ mesons is about the same 
for all reaction systems from C+C to Au+Au collisions 
as well as for all collision centralities and hence for very 
different densities inside the collision zone. The KN po- 
tentials on the other hand arepredicted to have a strong 
dependence on the density [Ty] . Yet, the data do not 
show such a dependence. This is a hint, that other effects 
besides the KN potentials might as well be important for 
the explanation of the different inverse slope parameters 
of the K + and of the K~ mesons. 

Comparisons of data from the KaoS-Collaboration to 
various transport-model calculations have been shown in 
several publications. Some of the more recent compar- 
isons can e.g. be found in [H, d, [45|, |4(| |47| and in the 
references therein. Another comprehensive report is in 
preparation [48[ . A comparison of the various transport 
models can be found in [491 ]. Here, we concentrate on the 



description of spectra and angular distribution and their 
sensitivity to in-medium modifications of kaons in dense 
nuclear matter. This comparison will be short as fur- 
ther developments are expected from the theoretical side. 
Most transport models use parameterizations for the KN- 
potentials that result in density-dependent K + and K~ 
selfenergies (see |l6[ for examples). A more recent con- 
cept to describe the in-medium properties via spectral 
functions of the K~ by Lutz et al. |l7| has not been 
implemented into transport models so far. Another ap- 
proach using coupled-channel G-matrix calculations [181 ] 
has been used for the K~N interaction in the Hadron 
String Dynamics Model (HSD) 0. 

In the following we compare energy distributions and 
polar angle distributions to results of calculations ob- 
tained with the Isospin Quantum Molecular Dynamics 
model (IQMD) 0, M§ and with the HSD model g|. 
These calculations have been performed by the authors 
of the respective codes analyzing the results within the 
experimental acceptance of the KaoS measurements. 

Figure |2"T1 displays the invariant cross sections for K + 
and for K _ mesons as a function of their kinetic energy 
for Au+Au collisions at 1.5 A GeV. In the upper part 
the comparison to IQMD calculations is shown, in the 
lower part the comparison to HSD results. In both cal- 
culations the energy distributions of the K + are slightly 
steeper than for the measured data. The yields favor the 
option with a repulsive K + N potential. For the K~, the 
slopes of both model calculations without an in-medium 
K~N interaction agree rather well with the data. The 
comparison of the absolute yields, however, does not al- 
low to draw a conclusion as the two calculations differ 
strongly. The discrepancy is only seen for K _ , while the 
results of the two transport models agree rather well with 
each other for the K + . 

As has been demonstrated, for the experimental data 
the inverse slope parameters of the K + are significantly 
larger than those of the K~. This trend is seen as well in 
the transport model calculations, rather pronounced in 
those calculations including in-medium KN interactions 
but already visible without them. 

Figure [22] compares the angular emission patterns of 
K + and of K~ in Au+Au at 1.5 A GeV, normalized to the 
yield at 8 c . m . = 90°, to the results of the transport models 
IQMD and HSD. For the K~ mesons from near-central 
collisions (0% - 18.1% cr r ) a rather flat distribution is 
observed, whereas the K + are preferentially emitted to 
forward and backward angles. The measured data are 
rather well described by both models. However, as can 
be seen in Fig. [2U this observable is hardly sensitive to 
the choice of the KN in-medium interaction. 

The differences in the inverse slope parameters and in 
the polar angle distributions of the K + and of the K~ 
mesons may as well be influenced by different emission 
times of the two particle species. This is demonstrated 
in Fig. |2"51 showing the results of IQMD calculations. The 
upper panel shows the density p reached in the collision 
zone normalized to normal nuclear matter density po as 
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FIG. 21: (Color online) Comparison of the energy distribu- 
tions of the invariant cross sections of K + and of K~ in inclu- 
sive Au+Au at 1.5 A GeV with two different trans por t- model 
calculations, IQMD [|] (upper part) and HSD (lower 
part). 
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FIG. 22: (Color online) Polar angle distributions of K + 
and of K~ in Au+Au at 1.5 A GeV. The lines are the re- 
sults of transport-model calculations, in the upper part using 
IQMD in the lower part using HSD 



high-density phase. According to IQMD the primary pro- 
duction of both, K + and K~ mesons, is isotropic. As will 
be discussed in detail in [48| the polar angle anisotropy 
of the K + is mainly caused by rescattering. The K~ on 
the contrary are emitted at a late stage of the reaction 
when the spectator matter has moved away and can not 
cause a significant anisotropy. 



a function of time. The lower panel shows the rate of 
emitted K + and K~ mesons as a function of their cre- 
ation time. It can clearly be seen that according to this 
transport-model calculation those K" leaving the reac- 
tion zone are created at a rather late stage of the reac- 
tion, significantly later than the K + which are mainly 
created during the high-density phase. This difference in 
emission times is caused by the strangeness-exchange re- 
action reabsorbing most of the K~ produced during the 



C. The Nuclear Equation of State 

Positively charged kaons are expected to carry informa- 
tion on the high density-phase of the heavy-ion collision 
due to two reasons: 

i) The K + production occurs via multiple collisions 
e.g. via processes like NN ^ NA followed by a subsequent 
interaction of the A resonance like NA ^ K + AN. These 
multi-step processes occur predominantly at higher den- 
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FIG. 23: (Color online) Results of IQMD calculations on the 
time evolution of central (6 = fm) Au+Au collisions at 
1.5 A GeV. Upper panel: density normalized to normal nu- 
clear matter density p/po as a function of the time. Lower 
panel: rate of emitted K + and K~ mesons as a function of 
their production time. 



sities Uli. 

ii) The rather large mean free path of the K + (ss 5 fm 
at normal nuclear density po) reduces the probability of 
further inelastic interactions prior to their emission. As 
a result of the K + production mechanism their yields are 
sensitive to the density reached in the collision which is 
related to the stiffness of the nuclear equation of state 
(EoS) as parameterized by the compression modulus 
defined as 



d 2 E/A(p,T) 



(dp) 2 



(7) 



which quantifies the curvature of E/A(p,T) at normal 
nuclear density p a . 

Figure [TT1 in Sect. lITTAl summarizes the multiplicities 
of K + mesons as determined in inclusive reactions of 
Au+Au, Ni+Ni, and C+C and normalized to the mass 
number A of the respective colliding nuclei as a func- 
tion of the beam energy. The K + excitation functions 
for all three collision systems rise strongly as expected 
due to the proximity of the threshold in NN collisions 
(Ethr = 1-58 GeV). The multiplicities per A are the 
higher the heavier the collision system is. This reflects 
that the K + mesons are predominantly produced in mul- 
tiple collisions which are more likely to occur the higher 
the density is in the reaction. 

Early transport-model calculations predicted that the 
K + yield in Au+Au collisions would be enhanced by a 
factor of about 2 if a soft rather than a hard nuclear 
equation of state is assumed [f| [5l| . Recent calculations 



take into account modifications of the kaon properties 
in the dense nuclear medium leading to a repulsive K + N 
potential which depends on the baryonic density [l6[ and 
which leads to a reduction of the calculated K + yields. 
To disentangle these two competing effects we use the 
ratio of two K + excitation functions [7|, one from C+C 
and one from Au+Au. The maximum baryonic density 
reached in Au+Au reactions is about 2-3 times nor- 
mal nuclear matter density while the increase in density 
in C+C collisions is significantly less pronounced. More- 
over, the maximum baryonic density reached in Au+Au 
reactions depends on the compression modulus of nuclear 
matter |5lll52T] whereas in C+C collisions this depen- 
dence is rather weak Q . Hence, the ratio of the K + mul- 
tiplicity per nucleon M/A in Au+Au to the one in C+C 
is expected to be sensitive to the compression modulus 
i^N- Furthermore, it provides the advantage that sys- 
tematic uncertainties within the experimental data are 
partly cancelled. This ratio in addition has turned out 
to be hardly affected by input quantities of the transport- 
model calculations which are less well known, like cross 
sections of individual reaction channels, the strength of 
the KN potentials, or the lifetime of the A resonance, as 
has been systematically studied in [9(. 

The upper and the middle panels of Fig. [24] show a 
detailed comparison of the K + excitation functions (a /A 
as a function of fibeam) for Au+Au and for C+C col- 
lisions with transport-model calculations. On the left 
hand side of the figure the data are compared to results 
from RQMD calculations [8j. On the right hand side a 
comparison to results from IQMD 0] calculations is pre- 
sented. The upper panels show calculations without KN 
potentials, the middle panels those with KN potentials. 
Solid lines denote calculations with a soft nuclear equa- 
tion of state (i^N = 200 MeV), dashed lines denote a 
hard EoS (K N = 380 MeV). 

The lower panels of Fig. [24] show the double ratio 
[M/A(Au+Au)]/[M/A(C+C)] as a function of the beam 
energy. Since due to the different energy loss in the Au 
and in the C target the effective energies for the K + 
production are slightly different, the fits as displayed in 
Fig. [11] were used for interpolation. The error bars con- 
tain the statistical uncertainties as well as those system- 
atic errors that do not cancel by calculating the dou- 
ble ratio (approximately 6%), added quadratically. The 
double ratios as determined from the various transport- 
model calculations are shown as well. Only the calcula- 
tions using a soft EoS are in agreement with the data. 

A reliable error estimate for the compression modulus 
i^N strongly depends on the transport model calculations 
and their input. A detailed study of this topic will be sub- 
ject of a theory publication [48j . The sensitivity of the 
double ratio [M/A(Au+Au)]/[M/A(C+C)} to the stiff- 
ness of the EoS within the standard version of IQMD is 
shown in Fig. [25] It compares the measured double ratio 
at 0.8 and 1.0 A GeV (shown as red bands) to results of 
IQMD calculations as a function of K^, both with and 
without KN potentials. The figure allows to establish up- 
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FIG. 24: (Color online) Upper and middle panels: Compari- 
son of the K + excitation functions (<r(K + ) /A) for Au+Au (red 
circles) and for C+C collisions (green squares) with RQMD 8] 
(left hand side) and with IQMD calculations [|| (right hand 
side). Solid lines depict a hard EoS, dashed lines depict a 
soft EoS. The upper panels show calculations without KN po- 
tentials, the middle panels calculations with KN potentials. 
Lower panels: Double ratio of the K + multiplicities per mass 
number M/A in Au+Au divided by the one in C+C and the 
comparison to the various transport-model calculations. 



per limits for K N , being < 180 MeV (at 0.8 A GeV) and 
< 240 MeV (at 1.0 A GeV) assuming a KN potential. 
Higher limits are obtained without KN potential result- 
ing in < 240 MeV (at 0.8 A GeV) and < 315 MeV (at 
1.0 A GeV). 

In Ref. [13] a detailed comparison of the K + produc- 
tion in the mass-asymmetric collision system C+Au as 
well as in Ni+Ni, both measured at the Kaon Spectrom- 
eter, with RQMD calculations is presented. Although the 
mean number of nucleons participating in the reaction is 
rather similar in both cases, the densities reached are 
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FIG. 25: (Color online) The double ratio 

[M/A(Au+Avl)]/\M/A(C+C)] calculated within the IQMD 
model (with and without KN potential) as a function of 
_Kn at 0.8 and at 1.0 A GeV. The experimental values are 
given as bands and allow to estimate upper limits for A"n as 
described in the text. 



significantly different. The compression in C+Au hardly 
exceeds the values obtained in C+C and therefore the 
K + yield in the calculations does not show a dependence 
on the stiffness of the EoS. For Ni+Ni the compression 
is significantly higher and the difference in the calculated 
yields is about 25%. Again, the RQMD calculation for 
A'n = 200 MeV is in good agreement with the data. 

In Ref. 0] the centrality dependence of the K + pro- 
duction has been suggested as a further observable to 
extract the stiffness of the nuclear EoS from heavy-ion 
collisions. As shown in Fig. QIJ] the K + multiplicity M in 
Au+Au at 1.5 A GeV rises stronger than linear with the 
number of participating nucleons ^4 pa rt- This is due to 
the production via multiple collisions which is needed to 
accumulate the necessary energy and due to the depen- 
dence of the number of collisions on the density reached 
in the reaction zone. In Sect. IIII Bl this rise has been 
quantified by a fit to the data of the form M ~ Ap art 
with a = 1.34 + 0.16. Figure compares this value 
(shaded area) to results of IQMD transport-model calcu- 
lations on the dependence of a on the compression mod- 
ulus AT N HHl. Only values for K N less than 250 MeV 
are compatible with the measured data assuming a KN 
potential, and K]$ < 320 MeV for the case without a KN 
potential. 
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FIG. 26: (Color online) Comparison of the exponent a. from 
the fit of M ~ Ap art to the K + multiplicities in Au+Au at 
1.5 A GeV (shaded area) with IQMD calculations as a func- 
tion of the compression modulus Ks ||- 



We would like to stress that two independent observ- 
ables, the centrality dependence of the K + yields, shown 
in Figs. [T51 and |2"61 as well as the system-size dependence 
presented in Figs. [TT1 and 1241 yield very similar results on 
the compressibility of nuclear matter when compared to 
transport-model calculations. Both observables support 
a soft nuclear EoS within the density regime explored 
by heavy-ion reactions at beam energies between 0.6 and 
2.0 A GeV. 



V. SUMMARY 

In this paper we have presented a comprehensive re- 
view of the production of charged kaons in heavy-ion col- 
lisions at incident energies (0.6 to 2.0 A GeV) below and 
at the respective thresholds in NN collisions as measured 
with the Kaon Spectrometer KaoS at GSI. This subject 
has been systematically studied by analyzing total pro- 
duction cross sections, energy distributions and polar an- 
gle distributions as a function of the size of the collision 
system, of the incident energy, and of the collision cen- 
trality. The key observations and trends can be summa- 
rized as follows: 

• The multiplicities of both, K + and K~ mesons, per 
mass number A of the collision system are higher in 
heavy collision systems than in light systems. This 
difference increases with decreasing beam energy. 

• The multiplicities per number of participating nu- 
cleons Ap ar t of K + and of K~ mesons within the 
same collision system rise stronger than linearly 
with j4 part whereas the pion multiplicity is propor- 
tional to Apart- Moreover the rise is rather similar 



for K + and for K - although the respective NN- 
thrcsholds for their production are significantly dif- 
ferent. 

• The K~/K + ratio is almost constant as a function 
of the collision centrality. At 1.5 A GeV this ratio 
is the same for Au+Au and for Ni+Ni collisions. 

• The inverse slope parameters of the energy distribu- 
tions of K + and of K~ mesons are higher in heavy 
than in light collision systems. 

• The inverse slope parameters of the energy distribu- 
tions of K + mesons are about 15 to 25 MeV higher 
than those of the K~ distributions. This is observed 
for all collision systems and for all centralities. 

• The polar angle distributions exhibit a forward- 
backward rise which is more pronounced for K + 
than for K~ mesons. K~ mesons produced in cen- 
tral collisions are emitted almost isotropically. 

From the systematics of these experimental results and 
from detailed comparisons to transport model calcula- 
tions the following conclusions on the properties of dense 
nuclear matter as created in heavy-ion collisions and on 
the production mechanisms of K + and of K - mesons can 
be drawn: 

• The K~ and the K + yields are coupled by 
strangeness exchange: Despite their significantly 
different thresholds in binary NN collisions the mul- 
tiplicities of K + and of mesons show the same 
dependence on the collision centrality. They are 
even similar for different collision systems. This 
can be explained by the K - being predominantly 
produced via strangeness exchange from hyperons 
which on the other hand are created together with 
the K + mesons. Strangeness exchange is predicted 
to be the main contribution to the K~ production 
in heavy-ion collisions at SIS energies by transport 
model calculations as well. 

• K + and K~ mesons exhibit different freeze-out con- 
ditions: Transport-model calculations predict dif- 
ferent emission times for K + and for K _ mesons 
as a consequence of the strangeness-exchange reac- 
tion. The K - are continuously produced and re- 
absorbed and finally leave the reaction zone much 
later than the K + mesons. This and the kinematics 
of the strangeness-exchange process manifests itself 
in an isotropic emission of the K - in central colli- 
sions and by systematically lower inverse slope pa- 
rameters of the K - energy distributions compared 
to K+ . 

• The nuclear equation of state is soft: The increase 
of M(K + )/A with the size of the collision system A 
points towards a dependence of the K + production 
on the density reached in the collision. The ratio 
of the K + multiplicities in Au+Au and C+C as a 
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function of the incident energy allows to extract 
the compression modulus of nuclear matter by 
comparing to transport-model calculations. Only 
calculations using a soft nuclear equation-of-state 
(i^N ~ 200 MeV) can describe the data. This con- 
clusion is rather insensitive to the various input pa- 
rameters of such calculations. A soft nuclear equa- 
tion of state is further supported by comparing the 
centrality dependence of the K + multiplicities in 
Au+Au collisions to transport model calculations. 

Our results demonstrate the importance of the 
strangeness-exchange reaction for the production and 
propagation of negatively charged kaons in heavy-ion col- 
lisions at incident energies from 0.6 to 2 A GeV, on the 
one hand coupling their yield to the K + -production, 
on the other hand causing a rather late emission of the 
K~. The production of positively charged kaons itself is 



strongly linked to the high-density phase of a heavy-ion 
collision allowing to determine that the equation-of-state 
of nuclear matter is soft within the density regime ex- 
plored by heavy-ion collisions between 0.6 and 2.0 A GeV. 
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